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The renormalized one-loop theory is a coarse-grained theory of corrections to the self-consistent field theory 
(SCFT) of polymer liquids, and to the random phase approximation (RPA) theory of composition fluctuations. 
We present predictions of corrections to the RPA for the structure function 5'(A;) and to the random walk model 
of single-chain statics in binary homopolymer blends. We consider an apparent interaction parameter Xa that is 
defined by applying the RPA to the small k limit of S{k) . The predicted deviation of Xa from its long chain limit 
is proportional to 7V~^/^, where A'^ is chain length. This deviation is positive (i.e., destabilizing) for weakly 
non-ideal mixtures, with XaN < 1, but negative (stabilizing) near the critical point. The positive correction 
to Xa fo"" low values of XaN is a result of the fact that monomers in mixtures of shorter chains are slightly 
less strongly shielded from intermolecular contacts. The depression in Xa near the critical point is a result of 
long-wavelength composition fluctuations. The one-loop theory predicts a shift in the critical temperature of 
0(N~^^^), which is much greater than the predicted 0{N^^) width of the Ginzburg region. Chain dimensions 
deviate slightly from those of a random walk even in a one-component melt, and contract slightly with increasing 
Xe- Predictions for S{k) and single-chain properties are compared to published lattice Monte Carlo simulations. 



I. INTRODUCTION 

Equilibrium properties of polymer mixtures that involve 
correlations over distances greater than a few Angstroms 
are rather well described by a family of closely related self- 
consistent field theories. Among these are Flory-Huggins 
(FH) theory for the free energy of a homogeneous mixture 
and the random phase approximation (RPA) theory of com- 
position fluctuations. Both of these are special cases of a 
more general self-consistent field theory (SCFT) of inho- 
mogenous polymer liquids. 

Theoretical treatments of polymer mixtures that go be- 
yond SCFT may be roughly classified as either microscopic 
liquid state theories or coarse-grained theories of correla- 
tion effects. Liquid state theories such as the polymer ref- 
erence interaction site model (PRISM)' and Born-Green- 
Yvon theories^ are designed to treat models with strong 
repulsive interactions, and attempt to predict all details of 
liquid structure from first principles. Coarse-grained theo- 
ries are based on more idealized models of chain conforma- 
tions and interactions, e.g., as continuous Gaussian chains 
with point-like interactions, similar to those used in coarse- 
grained SCFT. 

Over the past 20 years, a series of coarse-grained 
theories of correlation effects in polymer mixtures have 
made use of a functional integral representation of 
the partition function that was originally introduced by 
Edwards. ■'si^i^iiiSi^ All of these "field theoretic" studies of 
mixtures have thus far used an approximation in which the 
distribution of field fluctuations is approximated by a Gaus- 
sian functional. This approximation is referred to in the 
jargon of field theory (which we have adopted) as a "one- 
loop" approximation. An alternative way of deriving this 
theory, which we find more informative, is discussed in 
Sec. dVlb of this paper. 



A. Coarse-Graining and Renormalization 

The reason for studying correlations in coarse-grained 
models is the hope that, by doing so, one might de- 
scribe phenomena that are not captured by SCFT, but that 
are, in some sense, independent of short-wavelength de- 
tails. The first several studies of coarse-grained models of 
^jgjjjg3,4j5j6^7j8 y^gj-g^ however, plagued by a sensitivity of all 
predictions to the treatment of very short wavelength corre- 
lations, at length scales which the coarse-grained Hamilto- 
nian was never intended to describe. The way around this 
problem, which was first applied in this context by Wang^, 
is a renormalization procedure in which the sensitivity to 
monomer scale structure is absorbed into the values of a 
few phenomenological parameters. Wang"^ considered the 
one-loop prediction for the free energy of a homogeneous 
polymer blend, and showed that the sensitivity of this quan- 
tity to local fluid structure could be absorbed into an appro- 
priate renormalization of the Flory-Huggins interaction pa- 
rameter. More recently, waiS showed how a generalization 
of this procedure can be used to obtain unambiguous results 
for one-loop corrections to both collective and single-chain 
correlation functions at nonzero wavenumbers, using the 
same formalism for polymer blends and diblock copolymer 
melts. 

Development of this renormalization procedure required 
two insights that were absent from the interpretation of ear- 
lier studies of the one-loop theory: 

The first was that the SCF approximation should be iden- 
tified with the infinite chain limit of the one-loop theory, 
rather than a simple random-mixing approximation. Prior 
to Wang's work, the SCF approximation had often been 
conceptually identified with a simple mean-field approxi- 
mation that is obtained from a saddle-point approximation 
for the underlying functional integral. This mean-field the- 
ory ignores aU correlations except those arising from in- 



tramolecular connectivity, which is never a quantitatively 
useful approximation. Defining SCF parameters by a pro- 
cess of extrapolation to the infinite chain limit instead takes 
advantage of the fact that local correlations in liquids of 
long polymers depend only weakly on chain length, even 
in strongly correlated liquids. The difference between the 
free energy of a system of finite chains and that of a hy- 
pothetical reference system of infinite chains is largely the 
result of slight differences in long-wavelength correlations. 
These differences may be adequately described by a coarse- 
grained model. By treating the interaction free energy of 
an infinite-chain reference system as an input to the theory 
(which must be inferred from experimental or simulation 
data), one can obtain a renormalized theory that makes un- 
ambiguous predictions of corrections to SCF phenomenol- 
ogy- 

A second necessary insight was that the one-loop pre- 
diction for the structure factor S{k) at nonzero wavenum- 
ber k is affected by a renormalization of statistical segment 
length, as well as of the interaction free energy. The value 
of the statistical segment length 5 of a flexible polymer in 
a dense mixture is generally different than the "bare" value 
required as an input to the theory. The one-loop predic- 
tions for the magnitude of this difference**, are sensitive to 
short-wavelength correlations. This renormalization of b 
has no effect upon the predicted value of structure factor 
S{k) at exactly k — 0, and so did not need to be taken 
into account in Wang's interpretation of one-loop predic- 
tions for iS'(fc = 0). It becomes an unavoidable complica- 
tion, however, in calculations of single-chain correlations, 
and in applications to diblock copolymer melts, for which 
the dominant fluctuations occur at nonzero k. 

In our previous paper'", we focused primarily on clarify- 
ing the physical reasoning underlying our renormalization 
procedure, and on demonstrating that it yields a physically 
consistent interpretation. In this paper, we present predic- 
tions for binary homopolymer blends. Predictions for di- 
block copolymer melts will be presented elsewhere. 



B. Related Prior Work 

The analysis given here has several particularly close an- 
tecedents. 

Our analysis of the free-energy density of a homoge- 
neous blend and the long wavelength limit of S{k = 0), 
which is given in Sees. HJandlnll is closely related to that 
of Wang.'' The main differences arise from mathematical 
approximations that Wang introduced in order to obtain a 
completely analytic approximation for behavior near the 
spinodal. We have avoided all approximations other than 
the one-loop approximation itself, at the cost of some addi- 
tional numerical analysis. The additional rigor has allowed 
us to calculate some additional quantities (e.g., the absolute 
shift in the critical point) and to identify physical phenom- 
ena that are important far from the spinodal. 



Our analysis of corrections to Gaussian chain statistics in 
blends, given in Sec. |IV] is closely related to the recent the- 
oretical analysis of chain statistics in one-component melts 
by Beckrich, Johner, Semenov, and coworkers'' '^. Like 
us, but prior to our work, these authors introduced a renor- 
malized statistical segment length, defined by extrapolation 
to = oo, as part of a systematic calculation of small cor- 
rections to random-walk statistics. We initially failed to 
notice the connection to our own workr^ or to cite their 
work. These authors presented a beautiful analysis of uni- 
versal corrections to Gaussian statistics in an incompress- 
ible melt at intermediate wavelengths, less than the coil size 
but greater than a monomer size. They also obtained an- 
alytic results for one-loop corrections to the intramolec- 
ular correlation function at all wavenumbers in melts of 
equilibrium polymers, with an exponential distribution of 
chain lengths. By introducing a numerical procedure for 
extracting UV convergent parts of divergent integrals, we 
are also able to numerically evaluate one-loop predictions 
for melts and blends of monodisperse homopolymers and 
(elsewhere) for block copolymer melts. Here, we consider 
statistics of monodisperse polymers in binary blends. 

The most important antecedent of both our work and the 
recent work of the Strasbourg group is the much larger, and 
(generally) more sophisticated theoretical literature on cor- 
relations in semi-dilute polymer solutions. ''' At a technical 
level, the closest analog of our approach within this litera- 
ture seems to be the analysis of mono-disperse semi-dilute 
solutions by Ohta and Nakanishijiiii^ 

C. Outline 

The paper is organized as follows: Sec. II provides a 
brief review of the renormalized one-loop theory for the 
free energy density, and of our notation. Sec. Ill presents 
predictions for corrections to the RPA theory of composi- 
tion fluctuations. Sec. IV presents predictions for correc- 
tions to random walk statistics in polymer blends. Sec. V 
presents comparisons of theoretical predictions to results 
of published lattice Monte Carlo simulations. Sec. VI dis- 
cusses the physical relationship between the one-loop the- 
ory for the free energy and RPA theory of correlations. Sec. 
VII clarifies the physical origin of some features of the one- 
loop theory by relating them to corresponding features of 
the underlying RPA theory. 

II. THEORETICAL FORMULATION 

A. Model and Notation 

We consider a coarse-grained model of a binary ho- 
mopolymer blend. Let v be the average volume per 
monomer, for either species. Let the volume fraction, de- 
gree of polymerization, statistical segment length and pack- 
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ing length'* for species i be denoted by (pi, Ni, bi, and 
li = v/hf, respectively, for species i ~ 1 and 2. 

Let the potential energy be the sum of an intramolecular 
potential J/chain and a monomer-monomer pair interaction 



where 



f/int - 



\j dv j dr'C/,,(r-r')Q(r)c, (r') . (1) 



Here, Ci(r) = X]ms'^('' ~ R-miCs)) is the concentration 
of monomers of type i, in which Rmi(s) is the position of 
monomer s on molecule m of species i. The strength of the 
pair interaction is taken to be of the form 

f/.(r-r')^.(^:f^^ ^ + ^")5A(r-r') . (2) 

Here, (5a is a function (5a = F(A|r — r'DA^'^witha range 
A"'^ and an integral j dru\ = 1, so that Sa approaches a 
Dirac ^-function as A ^ oo. Here, i? is a dimensionless 
compression energy and xo is a dimensionless "bare" inter- 
action parameter The incompressible limit may obtained 
by considering the behavior when B is very large. 
Let 



(3) 



where Sci{r) — Ci(r) — (q), and let S'ij(k) = 
J dr e^^''^Sij{r) denote the Fourier transform. Let rii(k) 
be the intramolecular contribution to Su (k) from chains of 
type i. A continuous random walk model yields 



(4) 



where ^ — Nibf/6 is the radius of gyration for species 
i and g{x) — 2(e^^ — 1 + a;)/a;^ is the Debye function. 

Composition fluctuations in an effectively incompress- 
ible blend may be described by a scalar correlation function 



Sik) = Sn{k) = S22{k) ^ ~Si2{k) 



(5) 



Here, we follow Schweizer and Curro and our previous 
work, by expressing this scalar correlation function in 
terms of what we (and Wang'') refer to as an "apparent" 
interaction parameter Xa{k) , by defining 



s~\k) = n^\k) + n^\k) - 2vxa{k) 



(6) 



This is, of course, a generalization of the RPA expression 
for S{k), in which it is assumed that chains are random 
walks and that Xa{k) is independent of fc, N, and x^- 



B. One-loop Free Energy 



/id = -Tf- In 



N2V 



■ln( 



(8) 



is the ideal Flory-Huggins entropy of mixing, and 

Cmf = ^ J dr Uij{r)(ci){cj) 
= v~^Xo4'i4'2 



(9) 



is a simple mean-field approximation for the interaction en- 
ergy density. The remainder, denoted /con , is a correlation 
free energy. Here, and throughout the remainder of this pa- 
per except Sec. V, where we compare to simulation results, 
we set kT = 1. 

A one-loop approximation for /con has been obtained 
in previous work by considering a Gaussian approxima- 
tion for the distribution of fluctuations of a complex chem- 
ical potential field within the context of the Edwards auxil- 
iary field formalism. In Sec. IVII we discuss an alternative 
derivation in which this approximation is instead obtained 
by thermodynamic integration of the RPA correlation func- 
tion. However it is obtained, the resulting approximation 
for /coir is given, in the general case of a compressible liq- 
uid, by a Fourier integral 

/con = ^ [det|I + f2(q)U(q)|] (10) 

Here, the bold faced symbols fl{q) and U(k) indicate 2x2 
matrices with elements 6ijfti{q) and Uij{q), respectively, 
I denotes the identity, and det| • • • | is a determinant. Here 
and hereafter, we use the shorthand J^ = J d^q/ {2tt)^ for 
Fourier integrals. 

In the nearly-incompressible limit of interest here, Eq. 
[TOlfor /corr may be approximately reduced to 



/corr = ^ y" In [(ill + fl2 



2xovilin2) v] 



(11) 



To obtain the integrand of Eq. ( fTTT i. in a model with a pair 
interaction of the form given in Eq. we have both taken 
the limit _B 1 in Eq. ( fTOl ), and also replaced the Fourier 
transform of Sa by unity. This yields an expression for the 
integrand that is valid for |q| ^ A. In what follows, we 
regularize the one-loop approximation for /corr by treating 
A as a cutoff wavenumber, and simply restricting the inte- 
gral in Eq. (fTTI) to |q| < A. This is equivalent to the use of 
a pair potential in which the Fourier transform of (5A(r — r') 
is taken to be unity for q < A, and zero for q > A. The 
purpose of the renormalization procedure discussed below 
is to extract predictions that are independent of such details 
of the treatment of short-wavelength correlations. 



The Helmholtz free energy / per volume for this coarse- 
grained model may be expressed, in the incompressible 
limit, as a sum 



/ = /id + emf + /c, 



(7) 



C. Divergence and Renormalization 

It has long been known that, in a model of continuous 
Gaussian chains, the ideahzed limit A ^ 00 of point-like 
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interactions yields a one-loop expression for /con that di- 
verges as the cube of Ai^i^, In previous work, both Wang^ 
and wsiS have analyzed the structure of the high-^? diver- 
gence of this free energy density, and its physical signifi- 
cance. We found that Eq. (fTTI) for /corr can be expressed 



/corr — /corr ~^ f 



(12) 



of an "ultra-violet" (UV) divergent contribution /^^r' 
which increases with increasing A, and a contribution /* 
that is independent of A. The UV divergent contribution 
/^rr may be expressed as a sum'" 



A 



/bulk + /end 



(13) 



where 



.A 

/bulk 



.A 

J end 



1 



127r2 
3 



In 



121 

M 



2" 

+ 3 

M 

N2 



A^- 



A. 



6X0 hh 



A 



(14) 



and where I 



is proportional to l/iV, was identified'" as an excess free 
energy for chain ends. The quantity /i^^u^, which is inde- 
pendent of N, is instead a correction to the "bulk" interac- 
tion free energy density of a hypothetical reference system 
of infinitely long chains. 

To remove the explicit cutoff dependence, we re- 
interpret the one-loop free energy (excluding end effects) 



The contribution /i^j, which 



/ = /id + /i 



int 



/c^rd + r 



in which 



/ 



int 



Gmf + /bulk 



(15) 



(16) 



is identified as the interaction free energy of the relevant 
form of SCFT. The contribution is a non-universal 
chain end contribution that could (and generally should) be 
added to SCFT, to reflect the fact that the fluid structure is 
perturbed in the immediate vicinity of any chain end. The 
remaining cutoff-independent contribution /*, is identified 
as a universal correction to the SCF phenomenology. 

Eq. ( fT6b for /int simplifies considerably in the case of 
structurally symmetric polymers, with h = I2 = I- In 
this case /int is given, to within a composition-independent 
constant, by a function of the form 



/int = V ^Xe{A)(j>i(j)2 



where 



Xe(A) = 



1 - ^lA 



XO 



(17) 



(18) 



is an effective Flory-Huggins interaction parameter. This 
UV divergent expression for Xe was first obtained by 



Olvera de la Cruz et al? . The physical origin of the UV 
divergence of this quantity is discussed in Sec. IVIII 

We focus in the remainder of this paper on the UV con- 
vergent correction /*. It is convenient to introduce a no- 
tation for /* and related quantities, in which we use to 
denote a renormalized Fourier integral. The renormalized 
Fourier integral of an integrand A{q) is defined by a differ- 
ence 

f A{ci)^ I A(q)- f A{c^) (19) 

J q Jq J q 

between the unrenormalized integral A(q) and the UV 

divergent part of the integral, denoted by A{q). In 
a theory of continuous chains, the "UV divergent" part 
of Fourier integral is obtained by constructing a high-q 
asymptotic expansion of the integrand, in decreasing pow- 
ers of qR, and retaining only those terms in the expansion 
that lead to a UV divergent integralji^ 
With this notation, we may write 



1 

^ ^ 2 



In [(f^i + - 2vx^i^2) v] 



(20) 



The TV-dependence of /* may be isolated by non- 
dimensionalizing the renormalized integral in terms of a 
reference length R, which we take to be the root-mean- 
squared end-to-end R — Rgi of species 1 . It may be shown 
that all of the quantities in the integrand can then be ex- 
pressed as functions of qR and of a set of dimensionless 
parameters (f>i, xN, N2/N1, i?2/i?i. By this method, we 
find that /* may be expressed as a function of the form 



/* 



1 

R^ 



r{xN,4>i,N2/N,,R2/Ri) 



(21) 



in which /* is a dimensionless function given by the value 
of the non-dimensionalized renormalized integral. The re- 
quired set of dimensionless parameters is the same as that 
required by the standard non-dimensionalized SCFT for an 
incompressible blend. Note that, in Eqs. (|20] | and (ISTT i. we 
have written /corr as a product of a parameter x, but have 
not specified whether by x we mean xo. Xe. or (perhaps) 
the low-g limit of Xa- In fact, different choices lead to dif- 
ferent variants of the theory, as discussed in more detail 
below. 

Results for /* and related quantities that are reported 
here have all been obtained by numerically evaluating the 
relevant renormalized Fourier integrals. To do this, for each 
quantity of interest, we first numerically evaluate the un- 
renormalized integral with a large but finite cutoff |q| < A, 
and then subtract an analytic result for the UV divergent 
part of the regularized integral. The Fourier integrals re- 
quired by the one-loop theory may all be expressed as ei- 
ther one dimensional integrals with respect to a wavenum- 
ber q (for free energy density and its derivatives) or two- 
dimensional integrals with respect to q and a polar angle 
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(for correlation functions at nonzero wavenumber). Sub- 
traction of the UV divergent part from an unrenormalized 
integral yields an estimate of the renormalized integral that 
depends only weakly on the numerical cutoff A. Our final 
result for the value of a renormalized integral is obtained 
by repeating this procedure for several values of A ^ i?, 
and numerically extrapolating to A oo. 



D. Structure Function 

The long-wavelength limit of S{k) is related to the com- 
position dependence of the free energy density / by a sta- 
tistical mechanical theorem: 



where 



lim S ^(k) — v'^—TT . 
' " d<j>i 



fc^O 



(22) 



The quantity S ^ (0) is also related to the long-wavelength 
Umit of Xa(fc) by 



lim S-\k) ^ 2v{xs - Xa) 



(23) 



Here, and hereafter, we use Xa to denote the q ^ limit 
of the quantity Xa(k) defined in Eq. The quantity Xs 
is the spinodal value of Xa, given by 



2Xs 



1 



1 



(24) 



Using decomposition ( fTSl ) for /, and hereafter neglecting 
the explicit end effects, the one-loop approximation for Xa 
may be expressed as a sum 



Xa = Xe + X* 



in which 



Xe 



'2 d(l)l 



(25) 



(26) 



is the interaction parameter of phenomenological SCFT, 
and 



X 



(27) 



is a correction to SCFT arising from long-wavelength fluc- 
tuations. 

An explicit Fourier integral expression for x* may be 
obtained by differentiating the integrand of Eq. (|20] i for /* 
twice with respect to . The resulting approximation for 
X* may be expressed in terms of either Edwards' screened 
interaction '° or in terms of the correlation function 



S-\k-x)^^^\k) + n^\k)^2vx ■ 
When expressed in terms of S{k), it is found that 



(28) 



1 



1 1 



ry - l/02 + l/02_ (30) 

The quantities /i(0) and A(0) are proportional to the third 
and fourth derivatives of /id (01 ) with respect to 0i, respec- 
tively. Since ^{q), \{q) and rj are independent of x, the 
proximity to the spinodal in this expression is controlled 
by S{q) alone. 

By non-dimensionalizing the integral in Eq. ( |29] l, we 
may show that it yields a correction to A^Xa of the form 

Nx* ^^r{xN,<l>i,N2/N^,Rg2/Rgi) , (31) 

in which x* is a non-dimensional function that is defined 
by the non-dimensionalized integral. Here, 



N EE Nb^/v'^ 



(32) 



is an invariant degree of polymerization, in which (by con- 
vention) N — Ni and b — bi. 

In Eqs. ( |20] l. ( |2TI ). and ( [3TT i. we have intentionally ex- 
pressed /* and X* as functions of an interaction param- 
eter x^ without specifying whether this input parameter 
should be taken to be the bare parameter xo or some type 
of a renormalized value. In fact, different variants of the 
renormalized one-loop theory can be obtained by differ- 
ent choices for the interaction parameter used within these 
Fourier integrals. If we define /* by simply subtracting the 
UV divergent part of the one-loop expression for fcorr, we 
obtain a theory that is no longer UV divergent, but that re- 
tains an explicit dependence upon the unmeasurable "bare" 
parameter xo ■ Following the logic used to construct renor- 
malized expansions of field theories, however, we may also 
replace xo by a renormalized parameter within each inte- 
grand to obtain a theory from which all reference to the 
bare parameter have been removed. This substitution may 
be formally justified, within the context of a one-loop ap- 
proximation, by observing that the difference between ap- 
proximations obtained by using a bare or renormalized pa- 
rameter in the one-loop theory differ only at second order 
in a systematic loop expansion. 

In the remainder of this paper, we follow Wang^ by con- 
sidering a self-consistent one-loop approximation. In this 
approximation, also known as a "Hartree" approximation, 
X is replaced within the Fourier integral expression for x* 
by Xa = Xa(q = 0). In the resulting theory, Xa is given by 
a function of the form 



Xa = Xe +X*(01,Xa) 



(33) 



X'=-/ [^iHq)S^{q;x)-X{q)S{q;x)+v] (29) 



Here, Xe is related to the SCF (i.e., A^ oo) interaction 
free energy /int(0) by Eq. (|26] |. For purposes of numerical 



5 



evaluation, it is convenient to rewrite Eq. ( |33] l as an explicit 
expression Xe = Xa - X*('/'i,Xa) for Xe as an explicit 
function of Xa- 



III. COLLECTIVE FLUCTUATIONS 
A. Symmetric Blends, Critical Composition 

Figs. [Udisplays predictions of the renormalized Hartree 
theory for Xa in blends of structurally symmetric polymers 
(b ^ hi — 62 and N — Ni — N2) at their critical com- 
position (01 — 1/2). It shows predictions for XaN as a 
function of XeN blends with three different chain lengths, 
corresponding to TV = 64, 128, and 512. The difference be- 
tween XaN from XeN is largest for the shortest chains, and 
vanishes in the limit N ^ 00. The SCF predicts a critical 
point for a symmetric — 2. The actual critical point 
for such a blend is reached when XaN — 2. 

Fig. |2] shows the underlying dimensionless function 
X* = (XaN ~ XeN)VN as a function of XaN. Note that 
both Xa and %* vanish in the limit Xe ~ for any struc- 
turally symmetric mixture, with bi ~ 62 and iVi — N2, 
because this is an exact result for such a ideal mixture, in 
which the two species are physically indistinguishable. 



XaN 1 




FIG. 1: Self-consistently calculated versus XcJV for sym- 
metric binary blends with 0i = 1/2 for iV = 64, 128, 512 respec- 
tively (solid lines). The dotted line shows where Nxa = Nxe, 
which corresponds to the limit N ^ 00. Solid dots along the line 
XaN = 2 indicate critical points for finite chains. 



Notably, the dependence of x* upon XaN is not mono- 
tonic: For small values of XaN, the deviation x* = Xa~Xe 
is positive and increases with increasing Xa- Closer to the 
critical point, x* begins to decrease with increasing XaN, 
changes sign at XaN ~ 1.7, and is negative at the criti- 
cal point. Note that all the curves in Fig. ([U intersect at 
XaN = XeN ~ 1.7, where x* = 0. For any structurally 
symmetric blend, with bi = &2 and Ni = N2, we find that 



ife; -2 




XaN 



FIG. 2: The deviation A'^ ' N{xa — Xe) as a function of Xf^N for 
symmetric binary blends at the critical composition, <j)i = 1/2. 



the initial slope is 



d{x*N) {6/7T) 



lim 

XaN~.0 diXaN) 



3/2 



Nl/2 



(34) 



for any composition. 

At the critical point of such a blend ((j>i = 1/2 and 
XaN = 2), the one-loop prediction for x*^ approaches 
a negative limiting value 



lim x*N = -3.7N-^/^ 



(35) 



with a slope that diverges at the critical point. In the Hartree 
approximation, this yields a predicted critical value of 



(XeiV)c = 2 + 3.77V-i/2_ 



(36) 



If Xe is a decreasing function of temperature, this corre- 
sponds to a depression in the critical temperature Tc by an 
amount ST^ oc N^^/"^ 

The predicted decrease in x* near the critical point is not 
surprising: Near the critical point, long wavelength corre- 
lations decrease the level of local mixing between A and 
B monomers, thereby decreasing the apparent interaction 
parameter Xa and stabilizing the homogeneous phase. 

The origin of the positive values obtained for x* at small 
values of Xa is less obvious. In Sec. I VIII of this paper we 
show that this prediction is the result of a subtle depen- 
dence of local correlations in a polymer liquid upon overall 
chain length, which is discussed in much greater detail in 
an accompanying paper In a dense liquid of long polymers, 
the immediate vicinity of each monomer is crowded with 
monomers that belong to the same chain. This intramolec- 
ular "self-concentration" reduces the space available for 
monomers from other chains, thus creating a so-called cor- 
relation hole in the intermolecular correlation function. In 
a mixture, the correlation hole reduces the magnitude of 
the potential energy of interaction between monomers of 



6 



types 1 and 2. We show in the accompanying paper that 
the depth of the correlation hole, decreases slightly with 
decreasing chain length. For from the spinodal, this sim- 
ple packing effect causes Xa to increase with decreasing 
chain length, reflecting the fact that shorter chains are less 
strongly shielded from intermoecular contacts. Because we 
define Xe to be the limiting value of Xa in the limit of in- 
finitely long chains, this effect tends makes a positive con- 
tribution to the difference Xa{N) — Xe for any finite A'^. 
In the accompanying paper, we analyze this effect without 
reference to the one-loop theory, and derive the prefactor in 
Eq|34|by very different reasoning than that used here. 

A shift in the critical temperature by an amount STc oc 
iV~^/^ was predicted previously by Holyst and Vilgis^. 
The one-loop theory of Holyst and Vilgis did not use a sys- 
tematic renormalization procedure, but instead introduced 
an ad hoc cutoff A ^ 1/Rto avoid the UV divergence. 
This approach captured the correct scaling, but does not al- 
low a meaningful calculation of numerical values for one- 
loop corrections, because it yields corrections whose val- 
ues are sensitive to the exact numerical value chosen for an 
arbitrary cutoff length. 



B. Vicinity of the Spinodal 

We now discuss the behavior of the Hartree theory near 
the spinodal. This part of our analysis is similar to that 
given by Wang"^, who focused on behavior near the spin- 
odal, but who made several approximations that we avoid. 

A rigorous asymptotic expansion of the behavior of Xa 
near the spinodal may be obtained by expanding S{q;Xa), 
li{q) and X{q) in Eq. (|29] l around their values at g = and 
Xa^Xs- Let 



r 

T 



Xq) 

Xe)- 



The required asymptotic expansion of S{q) is 



(37) 



(38) 



where 



= Vl/{18cl>i<l>2lll2) 



By also expanding and X{q) around q 
an asymptotic expansion of the form 



(39) 

0, we obtain 



(40) 



where A, B, and C are dimensionless coefficients that are 
independent of r, but that that depend upon (pi, N2/N1, 
and R2/R1. Here, N is an arbitrary choice of a reference 
degree of polymerization, which could be taken to be A^i 
or N2 or an appropriate average of the two. 



The coefficients A and C in this expansion can be cal- 
culated analytically, by considering the singular behavor of 
the integral near g = as r ^ 0. The coefficient A is given 



1 



1 



327r \<j>iNi <j>^N2 



(41) 



Here, fi{q — 0) is the third derivative with respect to (pi 
of the ideal mixture free energy per monomer w/idl'/'i)- 
This quantity vanishes at the Flory-Huggins critical compo- 
sition, at which (piNi — 4>\N2, and so (as noted by Wang^) 
A{(j)i) also vanishes at the critical composition 0ic. At the 
critical composition, we find 



-1 



N 



(42) 



These values for A and C agree with corresponding results 
of Wang^. 

In the Hartree approximation, the value t^, of t (the re- 
duced "temperature") at the critical point is determined by 
the value of the constant B at the critical composition: The 
critical point occurs when r = at = 0ic, or when 



(43) 



If B is nonzero, an expansion of this form thus yields a 
shift in the critical value of XeN of order N^^^^. 

The coefficient B cannot be obtained from an asymp- 
totic analysis of the low-g behavior of the integral in the 
Eq.(l29]l. To understand why, recall that the problem ac- 
tually involves three length scales: the monomer size (or 
cutoff length A^^), the coil size R oc VN, and correla- 
tion length £, ^ £,q/ \fr. Very near the spinodal, £, ^ R. 
Our renormalization procedure removes all dependence on 
A^^, but leaves a dependence on both R and ^. Away from 
the critical composition, the integral in Eq. ( |29b develops 
an infra-red (IR) divergence of the form A/y/r as r — * 0. 
This divergence is the result of fluctuations with wavenum- 
bers q ~ and so its prefactor can be obtained by con- 
sidering the behavior of the integrand near q — 0. At the 
critical composition, the coefficient C can be isolated by 
considering the IR divergence of the derivative of Eq. (l29l l 
X* with respect to Xa, which then diverges as t^^^^. The 
actual value of x*, however, convergences even at the crit- 
ical point, and is dominated by fluctuations with qR ~ 1. 
As a result, the one-loop prediction for the shift in the crit- 
ical point can only be calculated by numerically evaluating 
the renormalized integral, as done here, without introduc- 
ing approximations for ili{q) and ft2{q) that are valid only 
for qR < 1. 

At the critical composition, the asymptotic expansion of 
the Hartree theory may be expressed more compactly as a 



Nr = N5t + 



Cv 



(44) 
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where 5t = t — t^, with C < 0. The theory exhibits 
strongly non-classical critical behavior for values of St less 
than a crossover value 6t*. Approximating St* by the 
value at which the y/r fluctuation correction is equal to St 
yields a reduced crossover temperature 

NSt* - TV^i , (45) 

as first noted by de Gennes.'^ As noted by Holyst and 
Vilgis^, the predicted 0{l/N) with of the critical region 
is thus much less than the predicted 0{N^^^^) magnitude 
of the shift of the critical temperature from its SCF value. 

Very close to the critical point, r ^ t*, the Hartree the- 
ory yields r cx t^, or S{q = 0) oc r"^. The predicted 
critical exponent of 7 = 2 is much larger than the known 
Ising critical exponent of 7 ~ 1.26. This is a well known 
defect of the Hartree theory. 

C. Off-Critical Blends 



1 



Xa/ Xb 



0.6 




^ 0.8^"-VX0.5 






y 1 


N = 4000 
1 1 1 1 



0.2 0.4 0.6 0.8 1 

Xe/ Xs 

FIG. 3: Self-consistently computed Xa versus Xe for off-critical 
blends with 61 = 62 and A^i — N2 and iV = 4000 at various 
values of (j) = (j>i. Both axes are reseated by the spinodal value 
Xs- Note that the spinodal line Xa/Xa = 1 is actually reached 
only at the critical composition of (pi = 1/2. The computed 
locus of turning points is depicted with dashed line. 



At any composition except the critical composition, the 
Hartree theory considered here does not have a spinodal: At 
off-critical compositions, because x* diverges like r~^/^ 
with decreasing r, there is no value of t for which r — 0. 
As noted by Wan^, in the Hartree approximation, this 
causes Xe to become a non-monotonic function of Xa or 
(equivalently) Xa to become a multi-valued function Xe- 
An example of this behavior is shown in Fig. [3] which 
shows predictions of renormalized Hartree theory for Xa 
vs. Xe in a series of off-critical blends of various composi- 
tions with 61 = 62, ^1 = N2, and N = 4000. The turning 
point at which Xe reaches its maximum value, which is also 
where S{q — 0) is maximum, is what Wang referred to as a 
"pseudo-spinodal". We do not attach any physical signifi- 
cance to this turning point, except as a particularly obvious 
failure of the Hartree approximation. 



The one-loop correction to Wi(k) in a homopolymer blend 
is given b y'°i'^ 

Su;,ik) = -^J dq^f (k,-k,q,-q)G,,(q) . (47) 

Here, Gii{q) is the screened interaction between any two 
i monomers, which is given in an incompressible blend at 
wavelengths g <C A by 



1 - 2vXe^in2n-^ 



The function ^f"^ = tjj^'^'' (k, — k, q, — q) is given by 
V'f^ = t^f^lk, -k,q, -q) -cj,(k)wj(q), 



(48) 



IV. SINGLE CHAIN PROPERTIES 

The one-loop theoryiSiJl also predicts slight deviations 
from random-walk statistics for individual chains. The 
one-loop approximation for the intramolecular correlations 
function is a form of renormalized perturbation theory in 
which each pair of monomers on a chain interact via Ed- 
wards' screened interaction, and in which the effect of this 
interaction is taken into account to first order in powers of 
the strength of the screened interaction'^. 

We consider the single-chain correlation function oJiCk). 
This is given for a Gaussian homopolymer by 



(46) 



where Rg i = Nibf/6 is the radius of gyration of species 
i and g{x) = 2{e^^ — 1 + x)/x^ is the Debye function. 



where uj^'^'' is the 4-point intramolecular coiTelation func- 
tion for an ideal Gaussian chain''*. It is straightforward 
to verify that ^aji(k) vanishes in the limit k = 0, as 
required to retain consistency with the requirement that 

If the high-g behavior of the above integral is regular- 
ized introducing a sharp cutoff at g = A, and using Eq. 
(|48l l for all q < A, the resulting integral is found to in- 
crease linearly with A as A is increased. We"^ and Beckrich 
et flZ. " have shown, however, that this UV divergence can 
be absorbed into renormalization of the statistical segment 
length. After this UV divergent part is subtracted, the re- 
maining UV convergent correction, which we will denote 
by 6u!* (k) may be used to characterize the deviation from 
random-walk statistics. More precisely, Sui* (k) is the devi- 
ation of the intramolecular correlation function (k) from 
that of a Gaussian chain with a renormalized statistical seg- 
ment length boo,i characteristic of a dense one-component 



8 



liquid of infinitely long chains. The correction (5tj*(k) is 
found to be smaller than uJi(k) by a factor proportional to 
iV^^/^, and exhibits a nontrivial dependence on k. 

The fractional deviation of the polymer radius gyration 
is given by (appendixlAli 



lim 



9(5tJi(k) 

a(fc2) • 



(49) 



where R^o,i = iV^fe^ j/6 is the prediction for a Gaussian 
chain with statistical segment length &oo,i- This fractional 
deviation is proportional to N^^^^. 



A. Monodisperse Melt 

We first consider the case of a monodisperse one- 
component melt, which we obtain by setting (j)2 — 0. The 
one-loop approximation used here reduces in the case of 
one-component liquids to that investigated by Semenov et 
^^ 11,12.20,21 ^ The only difference in our treatment of this 
case is our use of a numerical integration and renormaliza- 
tion procedure to obtain accurate results for monodisperse 
polymers. 

FiglH shows the predicted correction to the single chain 
correlation function in a monodisperse melt. From the low 
k behavior (see appendix|A]i we find that the predicted frac- 
tional change of the radius of gyration is 



Rio 



1.42 

Vn 



(50) 



We thus predict a contraction of finite chains, relative to 
a hypothetical Gaussian chain with a statistical segment 
length boo, in agreement with the conclusions of Ref. [see 
Eq. (20)]. 




N6u}{q) 



9-^9,0 

FIG. 4: Predicted correction to the single chain correlation func- 
tion (Debye function). 



B. Bidisperse Melt 

We next consider the behavior of a trace concentration 
of a species of chains of length TV (component 1) in a 



matrix of chemically homologous shorter chains of length 
A^2 = OiN (component 2). This hmit is obtained by calcu- 
lating the limit ^ of 6uJi (q). To represent chemically 
similar species, we set li = I2 and Xe = 0. This case re- 
duces to that of a mono-disperse melt when a = 1, and 
to a dilute solution of polymer in oligomeric solvent when 
a < 1. 



-0.5 




-1.5 - 



0.4 0.6 

a = N2/N1 

FIG. 5: The fractional change of the radius of gyration of an iso- 
lated chain of length A*' in a matrix of short chains of length aN, 
reseated with \/iV. 

In the previous subsection, we noted that the size of a 
single chain in the pure melt decreases as compared to the 
unperturbed value. It is also well known that a single chain 
swells in the dilute solutions with chemically similar small 
molecule solvents, which create a good solvent. In Fig|5] 
we show one-loop predictions for the fractional change of 
the radius of gyration of a test chain of length iV in a matrix 
of chains of length aN . The change is indeed negative for 
a close to 1 and becomes positive for a close to 0, passing 
through zero at a ~ 0.17. 

The limit a ^ 1 corresponds to a polymer solution in a 
good small-molecule solvent. The curve in Fig|5]diverges 
as 0.31/ a near a — Q (appendix|A|l. In this limit, where the 
matrix chains are point-like, Eq. ( l48T l for the screened in- 
teraction between a pair of monomers on a longer chain re- 
duces to a constant Gn (5) — vj (aN). Recall that the defi- 
nition of a "monomer" in an incompressible coarse-grained 
theory is simply the length of a chain that occupies an arbi- 
trary reference volume v. If, to simplify notation, the ref- 
erence volume is taken to be equal to the excluded volume 
vaN of a solvent molecule (i.e., a matrix chain), by setting 
aN = 1, we obtain an interaction Gii(g) = v character- 
ized by an excluded volume parameter equal to the solvent 
volume. Predictions of the theory in this limit are thus iden- 
tical to those obtained in the first order perturbation theory 
for a single chain with a point-like interactions v6{r)}^ It 
is straightforward to show that this perturbative result is 
consistent with the prediction^! that the longer chain will 
undergo substantial expansion if N2 ^ a//Vi, signalling a 
breakdown of first order perturbation theory. 
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C. Binary Blends 



V. COMPARISON TO SIMULATION 



We now consider binary blends with x ^ 0. For simplic- 
ity, we focus on symmetric systems, with Ni — N2 — N 
and &i = 62 = b. Fig. |6] shows the one-loop predictions 
for the fractional change of the radius of gyration squared, 
rescaled by Vn, as function of XeN for various composi- 
tions. Increasing x always causes both species in a blend to 
contract. The dependence of Rgi upon x is strongest when 
species 1 is a tracer in a matrix of 2 ((pi ^ 1). There is no 
dependence of Rgi upon x in the opposite limit (/)2 ^ 1 
of essentially pure homopolymer 1. In the dilute limit, 
4>i 0, the one-loop theory predicts a deviation SR^i that 
is a strictly linear function of x^- Furthermore, it can be 
shown that the slope of the limiting line is -0.62, which is 
-2 times that of the coefficient to the 1/a divergence for the 
bidisperse melt in the limit of small-molecule good solvent 
discussed in the previous subsection (appendixlAli. 

It is interesting to note that at the critical composition, 
(f)i = 1/2, the deviation induced by increasing x always re- 
mains less than the chain-length dependent deviation from 
ideal chain statistics that is already present in the one- 
component melt (x = 0). Also note that the radius of gyra- 
tion is a smooth, nearly linear function of x even near the 
critical point. This reflects the fact that the long-wavelength 
composition fluctuations that appear near the critical point 
have characteristic wavelengths much larger than the radius 
of gyration, and so couple very weakly to the conforma- 
tions of individual chains. 
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FIG. 6: The changes in the radius of gyration of chains at various 
compositions plotted as functions of in symmetric binary 
blends. The chains become more contracted as Xe^ increases 
and the effects is more pronounced for minority component. In 
the limit of (p — 1, the change becomes independent of Xe^' In 
the limit of <j) — 0, the change is linear in x^N, as can be told 
from the functional form of d, Eq.l|47|(. 



Here, we compare the predictions of the one-loop the- 
ory to the results of lattice Monte-Carlo simulation. We 
compare to bond-fluctuation model (BFM) simulations of 
symmetric blends at 0i = 1/2 Binder, Deutsch, and 
Miiller, and to BFM simulations of single-chain properties 
by Beckrich et a/. " . The bond fluctuation model is a lat- 
tice model on a simple cubic lattice in which each monomer 
occupies a set of 8 sites at the vertices of a cube. In all of 
the simulations discussed here, the monomer concentration 
is 1/16 monomer per lattice site, which is half the maxi- 
mum possible concentration in this model. The statistical 
segment length extrapolated for the infinitely long chain 
is 3.244, which has been reported in'^. In simulations of 
polymer blends, monomers of types i and j on neighbor- 
ing sites interact with an energy Uij, with uu — U22 and 
U12 > uii. Binder and coworkers have carried out semi- 
grand canonical ensemble simulations using two different 
variants of the rules for what constitute a "neighbor": One 
in which each monomer can interact only with monomers 
located at the 6 closest allowed sites, and another in which 
each monomer can interact monomers with monomers at 
any of the 54 closest allowed sites. 

A. Determining Xe 

To compare predictions of a coarse-grained theory to 
simulation data, one must somehow establish a relation- 
ship between the interaction parameters used in the simula- 
tions and the phenomenological interaction parameter used 
in both Flory-Huggins theory the one-loop theory. In an 
accompanying paper, we propose a way to do this, which 
we use here. The analysis given there refines and justifies 
a form of "modified Flory-Huggins theory" that was origi- 
nally proposed by Miiller and Binder.^'* There, we consider 
a perturbation theory for a lattice model of the type de- 
scribed above, in which free energy of mixing is expanded 
in powers of a perturbation parameter a = U12 — un. We 
identify the SCF interaction free energy per monomer fint 
by considering the behavior of the free energy of mixing in 
the limit of infinitely long chains. We show there that /i„t 
is given to first order in a power series in a by a function 

fint{e,^i)^kTxeM^-<f>i) (51) 

in which 

Xe=Z°°[ui2-Uii]/kT . (52) 

Here, z°° is an "effective coordination number" that is ob- 
tained by evaluating a value z{N) of the average num- 
ber of inter-molecular nearest neighbors per monomer in 
a one-component reference liquid, with a = 0, and ex- 
trapolating the value to the limit of infinitely long chains. 
The values of z°° used here were reported by Miiller and 
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Binder: = 2.1 for the model with 54 interacting sites, 
and — 0.307 for the model with 6 interacting sites. A 
more complete discussion of the reasoning underlying this 
prescription is given elsewhere. 

Our use of a first order perturbation theory to estimate 
fint in the one-phase region of the phase diagram is justi- 
fied by the fact that the critical value of a decreases as 1 /N 
with increasing chain length. The fractional error in our es- 
timate of Xe arising from the use of a first order expansion 
is thus expected to be of order 1/A^ near the critical point. 
For ^ 1, this error is thus expected to be much smaller 
than the corrections to Xa predicted by the one loop theory, 
which are of order N~^^^. The resulting 0{1/N) errors 
in our estimate of the "bulk" interaction parameter Xe are 
expected to be comparable in importance to the 0{1/N) 
corrections that arise from end effects. 



B. Composition Fluctuations in Blends 

Fig. Elshows the results of Mtiller for NS-'^{0) /2 plot- 
ted vs. XeN for bond-fluctuation model simulations of 
chains of length N — 6422.. These simulation used a variant 
of the BFM in which interactions extend over 54 neighbor- 
ing sites, and in which uu — —Ui2- Values of 5^^(0) were 
extracted from the semi- grand equation of state. We have 
calculated the SCF effective Xe parameter for each simula- 
tion from Eq. ( |52| | using the value of z°° ~ 2.1 reported 
by Mueller and Binder^'' 

In the upper panel of Fig. |7] the solid line is the RPA 
prediction and the dashed line is the one-loop prediction 
obtained using the appropriate value of iV = 290 (The 
value was reported as 240, which were computed from the 
statistical segment length of = 64 chains). The lower 
shows a comparison of simulation results for the difference 
X*N — {xa — Xe)N multiplied by to the universal 
curve shown in Fig. |2] We see that the self-consistent one- 
loop theory gives an extremely accurate description of devi- 
ations from SCFT far from the critical point, for XaN < 1, 
but appears to overestimate the stabilizing effect of criti- 
cal fluctuations. For this system, the suppression of Xa by 
long-wavelength fluctuations near the critical point is al- 
most exactly cancelled by the enhancement of Xa at smaller 
values of Xe, giving a critical value of XeN that is closer to 
the SCFT value than to the value predicted by the one-loop 
theory. 

Fig. m compares the predictions of the one-loop the- 
ory for the ration /Tc{N) for different values of N to 
the results of simulations by Deutsch and Binder^*. Here, 
= az°° N / {2k) is the critical temperature that is pre- 
dicted by SCF if we use Eq. (|52] | for Xe- The critical tem- 
peratures reported by Deutsch and Binder were identified 
for each chain length by a finite-size scaling analysis of 
simulations with different simulation cell sizes. Simulation 
results are shown for variants of the BFM model in which 
the neighbor interactions extend over 6 neighbors (circles) 




XaN 

FIG. 7: Comparison of the results of Miiller's bond-fluctuation 
model lattice simulations A'^S'^^(0)/2 (circles) to corresponding 
one-loop 1-loop predictions (dashed line), for A'' = 64 and iV = 
290, for a half-filled lattice"^. The calculation of values of Xe for 
simulation results is discussed in the text. 



and over 54 neighbors (triangles). For both versions of the 
models, /Tc{N) extrapolates to unity as iV — > oo. This 
confirms that first order perturbation theory prediction for 
Tc is indeed asymptotically exact in the limit N ^ oo. 

For both versions of the BFM, however, the measured 
deviations of /Tc{N) are significantly smaller than 
those predicted by the one-loop theory. This quantitative 
discrepancy must be, in part, a result of the inadequacy of 
the one-loop theory as a theory of critical phenomenon, and 
the need for a renormalization group approach near the crit- 
ical point. This data suggests that, if it is indeed possible to 
express universal corrections to /Tc{N) that arise from 
critical fluctuations as a power series in 1/ Vn, the coef- 
ficient of any 0(Af^^/^) correction must be much smaller 
than predicted by the one-loop theory. 

It is important to also note that the slight deviations of 
/Tc{N) from unity in Fig. |8]are not the same for the 
two variants of the BFM, and appear to of opposite sign. 
This model-dependence cannot be so easily explained as 
a failure of the one-loop approximation near the critical 
point, since we expect that a more sophisticated treatment 



11 



of critical phenomena in a simple coarse-grained model 
would also predict a universal dependence of /Tc{N) 
on the dimensionless parameter N. The observed non- 
universality could, however, be the result of corrections 
arising from contribution to fint that are second order in a, 
which are neglected in our estimate of Xe, and/or the result 
of end effects. Both of these effects are neglected in our 
analysis, and both are expected to yield model dependent 
corrections to T^/T^ that are proportional to 1 /N . It ap- 
pears to us that the dominant contribution to the deviation 
/Tc{N) — 1 for this data may well be a nonuniversal 
correction proportional to 1/iV (which appears to provide a 
better description of the N dependence, as originally sug- 
gested by Deutsch and Binder) rather than a universal cor- 
rection proportional to 1 /Vn. 




FIG. 8: Chain length dependence of the shift in the critical 
temperature. Comparison between 1-loop theory prediction and 
Deutsch and Binder's lattice simulatior?'. Simulation results are 
shown for variants of the BFM in which interactions extend over 
6 nearest neighbors (circles) and 54 neighbors (triangles). 



C. Single-Chain Correlations 

In Fig|9] we compare one-loop predictions and simula- 
tion results of Beckrich et al.^^ for the deviation 5u){k) of 
the correlation function for a monodisperse melt from the 
Debye function. Simulation results for 5 different chain 
lengths were taken from Ref . " . As noted by Beckrich et 
al, the collapse of data for from different chain lengths 
when ^/N5u}{q) is plotted vs. kRgQ already shows that 

the deviation is proportional to 1/ \/n, as predicted by the 
one-loop theory. Simulations and predictions agree quite 
well, though there do appear to be some small systematic 
discrepancies in the high k regime. The agreement is some- 
what better than that obtained by the Strasbourg group us- 
ing a more approximate treatment of a monodisperse melt. 
A similar level of agreement was obtained previously by 
that group in comparisons of the one-loop theory to sim- 
ulations of polydisperse equilibrium polymers, for which 
they obtained analytic predictions for (5a;(g) i"-'^i^°i^^'^^ . 
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FIG. 9: One-loop correction to the single chain correlation func- 
tion (Debye function). The data points are gleaned from Fig. 6 
of Ref.". Different symbols are lattice simulation (p = 1/16, 
b ~ 3.24) results for chains with different lengths: + (iV = 2048), 
o (1024), o (512), V (256), A (128). Nearfe = the data points for 
various chain lengths overlap and can not be distinguished from 
the original figure, so are symbolized as 



Miiller has examined the dimension of homopolymer 
chains of various lengths in binary blends with varying 
compositions and interactions^*^, using Monte Carlo sim- 
ulations of the BFM with the same parameters as the one 
used to study composition fluctuations^^. In Fig[TO]we re- 
plotted times his results for the fractional deviation 
of the mean-squared end-to-end vector for component 
A over a range of values of (f) = (t>A for a fixed value 
aN/ksT = 0.32, corresponding to XeN = 0.672, for 
6 different chain lengths, adapted from Fig. 5b in^*^. Here 
SR^ represents the deviation of R^ in a mixture from the 
ideal Gaussian chain value calculated using the statistical 
segment length of a sufficiently long chain. The original 
results were presented as the difference between the value 
in the blend and that in a melt of chains of the same length. 
We have switched the normalization by assuming that the 
radius of gyration measured in the melt is smaller com- 
pared to the ideal value by the amount 1A2/Vn. The fact 
that results for different chain lengths collapse confirms 
by itself that these deviations are proportional to 1/ ^/N. 
The method of obtaining the theoretical curve in FigfTOlis 
sketched in appendix|B] The qualitative trend of the simu- 
lation results were relatively well described with the theory. 



VI. THERMODYNAMIC INTEGRATION OF THE RPA 

The one-loop approximation has most often been pre- 
sented within the context of the Edwards formalism, as a 
Gaussian approximation for the functional integral repre- 
sentation for the partition function. Here, we discuss an al- 
ternative derivation that makes no reference to the Edwards 
formalism, in which the one-loop free energy is obtained 
by thermodynamic integration with respect to the strength 
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FIG. 10: Composition dependence of the fractional change of the 
mean squared end-to-end distance for homopolymer chains in a 
binary blends for systems with XeN = 0.672 and variable com- 
position. Simulation data were adapted from Miiller's data for 
Solid line is the one-loop prediction. 

of either the pair potential or the bare interaction param- 
eter, while using the RPA approximation for Sij{q). The 
derivation is closely analogous to the usual derivation of 
the Debye-Hiickel approximation for excess free energy of 
an electrolyte solutions, which is obtained by combining an 
expression for the free energy as an integral with respect to 
a "charging" parameter with the Debye-Hiickel linear self- 
consistent field approximation for the correlation function 
(see, e.g., McQuaiTie^'', pgs. 328-339). 



where 

5,fA(q) = [J7-l(q) + [/(q)]r.i (57) 

is the RPA (or mean-field) approximation for Sij (q) in a 
compressible liquid. We thereby find that the one-loop ap- 
proximation for the free energy density is thermodynami- 
cally consistent with an RPA approximation for Sij (q). 

It follows that the one-loop free energy may be obtained 
by a thermodynamic integration similar to that used to ob- 
tain the free energy in the Debye-Hiickel theory of elec- 
trolytes: Consider a sequence of systems with a reduced 
pair potential XUij, with a charging parameter A in the 
range < A < 1. The one-loop correlation energy (which 
vanishes in the A = ideal gas state) may be expressed as 
an integral 

/corr ^Ijy^l^ {q)S^/^ (q; A) (58) 

in which S^^^{q; A) is the RPA correlation function for a 
compressible mixture with a reduced pair potential XUij . It 
is straightforward to confirm that Eq. ( fTOl l is recovered by 
explicitly completing the integral with respect to A. 



B. Variation of xo 



A. Variation of the Pair Potential 

Our analysis starts with the following exact expression 
for the dependence of the free energy F = — In Z upon the 
pair potential Uij (r — r') in any model with a pair poten- 
tial: The perturbation SF in the free energy induced by an 
infinitesimal perturbation 5Uij (r — r') in the pair potential 
is given by an integral 



SF 



1 



dr / dr' (c,(r)cj(r')> SU,j{r - r') (53) 



The corresponding perturbation of the correlation free en- 
ergy density fcon may expressed as an integral 

Sfcorr^^l J dv S,,{V) 5U,j[v) (54) 

or as a Fourier integral 

<5/corr = \ j 'S'y (q) ^U^j{q) , (55) 

where 5„(r) = (c,(r)cj(0)) - (q)(cj). 

By explicitly evaluating the corresponding variation of 
Eq. ([Tol l for the one-loop correlation free energy for a com- 
pressible liquid, we find that, in this approximation, 

Sfcorr^l J S,fA(q)5C/.,(q) (56) 



For many purposes, it is useful to consider the simple 
special case of a symmetric mixture of structurally identi- 
cal chains, with bi = 62 and iVi = iV2, with a pair po- 
tential of the form given in Eq. (|2|i. In this case, the limit 
Xo = yields a true ideal mixture, for which /con is inde- 
pendent of composition. In a nearly-incompressible sym- 
metric mixture, with Bv ^ 1, the composition dependence 
of /con may thus be obtained by starting from a nearly in- 
compressible reference state with xo = 0, and calculating 
how fcorr varics as xo is increased. As a special case of 
Eq. ( fSST l, it is straightforward to show that 



S/c, 



dxo 



dr 5i2(r)5A(r) 



(59) 



If we take the incompressible limit, for which 5*12 (r) = 
— S'(r), and the limit A — > 00 of point-like interactions, 
Eq. i6l[ reduces to 



or 



lim — = ~v lim S{r) 

S/corr f or \ 

dxo 7q 



(60) 



(61) 



in Fourier space. 
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The one-loop theory may be obtained by using the in- 
compressible RPA approximation for S{q) in Eq. ( |6T] ). In 
a symmetric blend, the incompressible RPA yields 



Using this assumption, it is straightforward to show that 

Suiq) = -0i02f^('z) and that Sii{q) = +(l)i(f>2^{q) in 
any incompressible ideal mixture. 



S 



RPA 



(62) 



An idealized IRPA model of Gaussian chains with point- 
like interactions, however, yields a UV divergent integral 
in Eq. ( |6TT i. This may be regularized by restricting the 
integral to g < A. The regularized (but unrenormalized) 
one-loop free is then given, to within a composition inde- 
pendent constant, by an integral 



Jo Jq 

where the Fourier integral is constrained to |q| < A. 



(63) 



VII. CORRELATION HOLE EFFECTS 

Several features of the one-loop approximation for fcon 
may be understood by relating them to corresponding fea- 
tures of the incompressible RPA description of the so- 
called correlation hole. To show this, it is useful to con- 
sider the behavior of the derivative dfcorr/dxo in the sim- 
ple case of an ideal statistical mixture, which is obtained in 
our model by setting bi = b2, Ni ^ N2, and xo = 0. In 
this limit, Eq. ( |62] i for S{q) reduces to 



S{q) = M2^{q) 



(64) 



This expression can also be obtained more directly by the 
following physical arguments: Let 5(r) = {5c{r)6c{{))) be 
the overall correlation function in a one-component refer- 
ence liquid of either pure 1 or pure 2, and let H{r) and ri(r) 
be the inter-molecular and intramolecular correlation func- 
tions, so that S{r) = Vl{r) + H{r). An ideal mixture may 
be created by randomly labelling a fraction 0i of the chains 
in such a one-component liquid as species "1", and the re- 
mainder as species "2". In the resulting ideal mixture, the 
probability that a particular pair of monomers on different 



chains will be on chains labelled i and j is simply 1 



It 



follows that the inter-molecular correlation function Hij (r) 
in such a mixture is simply Hij{r) — (f)i(pjH{r). Sim- 
ilarly, because the probability that two monomers on the 
same chain belong to a chain of type i is (pi, the intra-chain 
correlation function ^ij{r) in such a mixture is ^ij{r) = 
6ij(piil{r). Thus far, everything we have said is exact. 

In the idealized limit of a nominally incompressible liq- 
uid, however, we assume that S{r) ^ 0, and thus that 
f2(r) + H{r) = 0. In an incompressible liquid, the inter- 
chain correlation function H{r) must thus have a "corre- 
lation hole", which is a mirror image of the intra-chain 
correlation function il{r). This implies that, in an in- 
compressible statistical mixture, Hij{q) = —(j)i(l)ji}{q). 



A. UV divergence (Unrenormalized Theory) 

The UV divergence of the correlation correction to Xe in 
Eq. JTSl l is the result of a pathological feature that arises 
from a combination of the assumption of rigorous incom- 
pressibility, at all wavelengths, in a model of continuous 
Gaussian chains. In the limit xo = and A 00, Eqs. 
( |6T| l and ( |64l l may be combined to show that 



dxo 



v4>i^2H{r^0)^-v^i^2n{r = 0) , (65) 



In the limit A 0, the predicted free energy is sensitive 
to the value of H{r) only at ?■ = because the range of 
the pair interaction Uij vanishes in this limit. It is well 
known, however, that the intramolecular correlation func- 
tion fl{r) for a continuous Gaussian thread diverges as 
il(r) cx 1/r as r ^ 0. Recall the origin of this divergence: 
The number of monomers within a region of volume 
near a test monomer is of order the number of monomers 
g in a segment of chain of size r — giving an aver- 
age concentration g/r^ ~ l/{rb^) within a region of size 
r, or ri(r) oc l/r. When combined with an assumption 
of rigorous incompressibility, however, the divergence in 
51(r) yields an infinitely deep correlation hole in H{r) at 
r = in the one-component reference liquid. Taken lit- 
erally, the incompressible RPA would yield negative val- 
ues of (c(r)c(O)) in a region around r = 0. This leads 
to the nonsensical prediction of a divergent reduction in 
H{r = 0), and thus to a divergent negative correction to 

Xe- 

The UV divergence obtained from the IRPA that under- 
lies the one-loop theory could be avoided by the use of any 
model that avoids the use of any of: i) a point-like inter- 
action, ii) a continuous Gaussian thread model, or iii) rig- 
orous incompressibility. For example, in the thread limit 
of the PRISM theory^", one retains a model of continuous 
Gaussian chains, and a point-like effective interaction (i.e., 
direct correlation function), but avoids complete nonsense 
by allowing for some compressibility at short wavelengths. 



B. A'-dependence (Renormalized Tlieory) 

In the renormalized one-loop theory, the UV divergent 
part of the correlation correction to Xa is absorbed into the 
definition of Xe(A). Once the divergence is removed, the 
renormalized theory can make concrete predictions about 
the difference x* = Xa — Xe- Because the one-loop ex- 
pression for X* is proportional to 7V^^/^, Xe may be inter- 
preted as the limit of Xa extrapolated to the limit N ^ 00. 
We now show that the value of x* is sensitive to the slight 
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difference between the depth of the correlation hole in a 
system of finite chains and that in a hypothetical system of 
infinite chains. 

To show this, consider the derivative of Eq. ( |20l i for 
f*{Xa) with respect to Xa, evaluated (for simplicity) in a 
reference state with Xa = Xo = 0. A straightforward dif- 
ferentiation yields 

df* 



1™ a 

Xa^O dXa 



-vM2 / niq) . (66) 

The only differences between Eq. ( l66b and the limit Xo = 
of Eq. i6T[ for 9/corr/<9xo are: i) the use of a renormal- 
ized integral (denoted by an asterisk) to remove the UV 
divergence, ii) the use of Xa rather than xo as the input 
parameter in a Hartree approximation. 

The UV divergent part of the integral in Eq. ( l66b is given 
by a difference 

df 



lim 

Xa-*0 dxa 



v'M2 / [niq)-n^iq)] 



(67) 



in which 0,-^{q) is an asymptotic approximation for i}{q) in 
the limit q ^ R^^. This quantity is given by 

n^{q) EE l2/{vq'^b^) . (68) 

Let rt{q; N) denote the intramolecular correlation function 
for a single-component melt of chains of length N. Be- 
cause il^{q) is an asymptotic expansion of il{q; N) in the 
limit qR ^ 1, it may be interpreted equally well as either: 

i) The high-g behavior of N) for chains of fixed 
length N, or 

ii) The asymptotic behavior of N) at fixed q in the 
limit N ^ oo. 

By adopting the latter interpretation, and inverting the 
Fourier transform, Eq. (|67| i, may thus be expressed as a 
difference 

df* 
lim - — 

Xa^O dXa 

or, equivalently, 

lim = t)2 lim \H(r; N) - Hlv; (x)] (70) 

Xa^O dXa 

where H{y\ N) denotes the value of the inter-molecular 
correlation function H{y) in a one-component melt of 
chains of length N . 

The renormalization procedure that was originally intro- 
duced remove the UV divergence is thus seen to be equiv- 
alent to the subtraction of the N ^ oo limit. The resulting 
difference relates x* — Xa — Xe to the difference in the 
depth of the correlation hole at r = in a system of finite 
chains and that in a hypothetical system of infinite chains. 
The renormalized coarse-grained theory correctly predicts 
the behavior for x* for small values of Xa because it cor- 
rectly describes this subtle difference between finite and 
infinite chains, despite the fact that it predicts a UV diver- 
gence H{r; N) itself. 



w0i(/i2 lini[r2(r;iV) - f7(r;oo)] (69) 



VIII. CONCLUSIONS 

We have use the renormalized one-loop theory to analyze 
corrections to the RPA description of composition fluctua- 
tions in binary blends, and to the random walk model of 
chain conformations. 

Our treatment of long-wavelength composition fluctua- 
tions is closely related to an earlier treatment by Wang, who 
focused primarily on behavior near the spinodal.'' Our pre- 
dictions can be summarized by considering the difference 
X* = Xa — Xe between an apparent interaction parameter 
Xa that would be inferred from the magnitude S{q — 0) 
of long-wavelength composition fluctuations and a param- 
eter Xe that we identify as the SCFT interaction parameter, 
which is the limit of Xa as oo. 

We show that the deviation from SCFT can be under- 
stood as the result of two competing physical effects: In 
the limit Xe^ <C 1 of a nearly ideal mixture, Xa is greater 
than Xe because of a simple packing effect: Monomers on a 
finite chain are slightly less strongly shielded from contact 
with monomers of other chains than are monomers on an 
infinite chain, and are therefore more exposed to unfavor- 
able interactions with chains of the other species in a binary 
blend. We discuss this effect from another point of view in 
an accompanying paper. Closer to the spinodal, however, 
the build up of long- wavelength correlations causes the dif- 
ference Xa — Xe to decrease with increasing Xe- Near the 
critical point, both effects yields contributions to Xa that 
are proportional to iV~^/^, and comparable in magnitude. 
An accurate prediction of the magnitude and sign of the 
absolute shift in the critical value of Xe^ thus requires an 
accurate description of both of these competing effects. A 
simple renormalized Hartree theory predicts an overall en- 
hancement in the critical value (A^Xe)c (or a depression in 
Tc) for a symmetric blend by an amount proportional to 
7V~^/^, as first suggested by Holyst and Vilgis,^'^ 

Comparison to the data from lattice Monte Carlo simula- 
tions of symmetric blends at the critical composition indi- 
cate that the one-loop theory gives an excellent description 
of the packing effect that dominates corrections to the RPA 
for Xe^ ^ 1, but a mediocre description of correlation 
effects closer to the critical point. It appears that the self- 
consistent one-loop (or Hartree) theory gives a reasonable 
semi-quantitative description of deviations from SCFT, but 
tends to overestimate the magnitude of the depression in Xa 
near the critical point, and thus to predict too high a value 

for (Xe^)e. 

Predictions for corrections to random-walk statistics 
were also compared to the available simulation data. Pre- 
dictions for deviations of the intramolecular correlation 
function from the Debye function for chains in a monodis- 
perse single-component melt agree very well with the sim- 
ulation results of Beckrich ef aZ. , " and deviate only slightly 
from a more approximate theoretical treatment of monodis- 
perse polymers by the same group. Predictions of the com- 
position dependence of the radius of gyration in a binary 
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blend with Xe 7^ show a similar level of quantitative 
agreement with the results of lattice Monte Carlo simula- 
tions. 



APPENDIX A: RADIUS OF GYRATION 

The fc-dependence of the correction to the intramolecu- 
lar correlation function of a homopolymer, Eq.(|47l). comes 
from that of t^I^^ (k, — k, q, — q). A straightforward Taylor 

expansion of tAI^^ shows that in the low k regime, to second 
order in an expansion in powers of |k|. 



where 6 is the angle between vectors k and q, and 

/(^) 



(Al) 



+ 



3 

1 



3 V a:' 
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x^ 


x^ J 



(A2) 



This expansion has been obtained previously in the studies 
of chain sizes in semi-dilute solutions'"* and in melts.-" Us- 
ing the fact that Lu{k) ~ Nf{l - PRl ^S) in the low-k 
region and Eq.(l47b. one finds: 



(A3) 



This is the general expression used to calculate the frac- 
tional change of the radius of gyration. 

Consider the variation of R^ for a tracer amount of poly- 
mer species 1 in a matrix of chains of species 2. In the rel- 
evant limit 01 0, the screened interaction Gu reduces 
to: 



1 



g{q^Rg^a. 



(A4) 



where g{x) = 2(e^^ — 1 + a;)/x^ is the Debye function 
and i?go,2 = N2i>2/Q- Substituting this expression into 
Eq.( IA3l ). and hereafter taking 61 = 62 = we find: 



-5R 



3,1 



F{a) 



p2 



a 



2F(0) • XeiVi 



where 



63/2 r- 



dy 



fiy')y' 



(A5) 



(A6) 



where a = N2/N1, and where y = q^R^Q i- F{a) is 
convergent when a = and is UV-divergent for any non- 
vanishing a. In the latter case, a constant 4/3 has to be 



subtracted from the integrand to account for the renormal- 
ization of the statistical segment length. 

Three special cases of this result are noted in the main 
text: The fractional shrinkage of R^ in a homopolymer 
melt, which is discussed in subsection II V Al is obtained 
by setting Xe = and a = 1. This yields vN SR^/Rgfl = 
F{1) = —1.42. In section HV B I the curve shown in Fig. 
|5] is given by F{a)/a. Its asymptotic behavior in the 
limit of a small molecule solvent, a ^ 0, is given by 
F{0)/a = 0.31/q!. In section ITVCl the Xe^^ dependence 
of the fractional change of Rg in the limit of is 

gwenhyd{y/NdRl/RlQ)/d{xeN) = -2F(0) = -0.62. 



APPENDIX B: MEAN-SQUARED END-TO-END 
DISTANCE 



We calculate the one-loop correction to the mean- 
squared end-to-end distance R^ using the method described 
in^^ (appendix III of chapter 5). There, the same quan- 
tity was calculated for a finite segment of an infinitely long 
chain in a semi-dilute solution, using a Lorentzian approx- 
imation for the screened interaction. The analogous result 
for a chain of species 1 in an incompressible binary blend 
is obtained by using the appropriate screened interaction 
Gij. The final result for a symmetric blend (A^i — N2 and 
bi = 62) is 



6Rl 
Rl 



24 r %2)y4 i^^^vg(y2) 

ay i _ Ar_/_9N (B1) 



5(2/2) 1 _ 2</)i02Xe^5(2/') 



where y = q^R^Q, and 



h{x) 



1 4 6 \ ^ 2 6 



(B2) 



This expression reduces to Eq. (15) of Ref4^ for a one- 
component melt (02 — 0). 

The integral is divergent in the high-g region, and the 
dangerous term in the integrand is "1". This is precisely 
the divergent term identified by Edwards, which he iden- 
tified with as a renormalization of the statistical segment 
length. We subtract this term from the integrand in order 
to obtain a UV convergent prediction for 6R^, which is the 
deviation of R^ from the value Rq (N) predicted by a Gaus- 
sian model with a statistical segment length equal to that 
obtained for infinitely long chains. For a pure homopoly- 
mer melt, by numerically completing the convergent in- 
tegral, we recovered the result reported in'^: SR^/R^ — 



-1.40/ViV. 
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